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SECTION  I 

INTRODUCTION 

This  thesis  contains  a  theoretical  investigation  into  non-linear 
electromagnetic- field  theory  as  applied  to  third  harmonic  generation 
in  bulk  materials.  Such  a  topic  is  of  interest  since  the  invention 
of  the  field  known  as  non-linear  optics.  With  the  high  power  levels 
available  from  lasers,  practical  experiments  in  non-linear  optics 
have  become  a  reality.  Of  ccurse,  the  generation  of  harmonics  and  the 
heterodyning  of  different  frequencies  or  wavelengths  is  desirable 
for  various  reasons.  The  most  compelling  reason  for  investigating 
third  harmonic  generation  (THG),  the  topic  of  this  thesis,  is  for 
use  as  a  tool  in  investigating  the  structure  of  ultra-short  laser 
pulses  known  as  pico- second  pulses. These  pulses  are  on  the  order  of 
pico-seconds,  a  time  regime  so  short  as  to  be  unresolvable  in  real 
time  by  present  techniques.  As  a  result,  non-linear  processes  such 
as  harmonic  generation  are  used  to  infer  the  characteristics  of  these 
pico-second  pulses.  It  is  for  this  purpese  that  the  investigation 
was  undertaken. 

In  format,  this  thesis  consists  of  introduction,  the  development 
of  the  THG  theory  and  several  examples  of  the  application  of  the 
theory  to  elementary  fundamental  functions  such  as  a  sine  wave,  a 
gaussian  pulse  and  a  linearly  chirped  pulse.  The  theory  begins  vdth 


1.  DeMaria,  A.  J.,  D.  A.  Stetser,  and  W.  H.  Glenn,  Jr.,  "Ultrashort 
Light  Pulses,"  Science,  Vol.  156,  No.  3782,  pp  1557-1568,  23  June 
(1967). 
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Maxwell's  Equation  and  concludes  with  an  expression  for  THG  under 
conditions  of  a  uniform  pxane  wave  in  an  isotropic  non-linear  medium. 
As  a  result  of  the  developed  theory,  the  influence  of  excitation  and 
material  parameters  can  be  evaluated.  Appendices  are  included  to 
elaborate  on  portions  of  the  various  developments  that  would  be 
unwieldy  in  the  body  of  the  thesis  as  well  as  to  elaborate  on 
interesting  issues  of  insufficient  importance  to  be  included  as  part 
of  the  principal  effort. 
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SECTION  II 

A  THEORY  OF  THG 


A.  Definitions.  This  section  of  the  paper  presents  a  descrip¬ 
tion  of  TUG  developed  for  a  uniform  plane  wave  in  on  isotropic 
non-linear  medium.  The  equations  that  are  obtained  are  sufficiently 
general  to  take  any  real  temporal  function,  whether  casual  or  not, 
and  obtain  the  resulting  third  harmonic  conversion  in  both  the  time  and 
frequency  domains.  For  a  point  at  which  to  begin,  the  definition  of 
a  Fourier  transform  pair  is  given  and  alao  symbology  that  will  be 
employed  in  the  mathematical  developments.  Thus,  the  Fourier  trans¬ 
form, 

09 

mCO 

and  its  inverse  transform, 


J  F((o)e'^i)tdw  , 

mUO 

are  given.  In  the  functional  notation  to  be  used,  functions  designated 
by  script  capital  letters  belong  in  the  time  domain  and  those  designa¬ 
ted  by  block  capital  letters  are  in  the  frequency  domain.  Also,  the 
Fourier  transform  relationship  may  be  symbolized  as 
3{"t)< — ►  F(oj) 

The  functions  of  interest  arc  defined  as  the  fundamental  function 
or  exciting  function,  £f(z,t),  as  , 

mCQ 


k 


the  generated  function  or  response  function,  £  (z,t),  os 

CD  S 


eg(*  ,o 


)e-Jk((o)zeMda)  # 


ond  the  polarizability,  (P(z,t),  as 
(P(z,t)=v-  I  P(z  >a,)e1u)tdij  . 


Ef(z,„)  ond  E^(z,«o)  represent  their  respective  temporal  functions 
leas  phase  factors  associated  with  propagation  while  P(z,«)  la  a  com¬ 
plete  description  at  all  points. 

B.  Relating  E  (z,»)  to  P(z.m).  Beginning  with  Maxwell's 

O 

Equations,  a  matliematical  development  will  be  pursued  thot  results  in 
finding  E  (z,«)  as  a  function  of  P(z,w).  There  will  be  approximations 

O 

made  os  the  development  progresses  and  these  will  be  explicitly  noted. 

Maxwell's  Eouotions: 

— x 


VXE  .-^(T 

^  *=  p 

V  -<B  «  o 


and  the  additional  relationships: 

^  n>r 

x*  — *  • 

J?”  €of  +<P 

Using  Maxwell' 8  Equations,  a  wove  equation  is  derived.  It  Is 
assumed  that  the  medium  is  charge  free  (P«0)  ond  of  zero  conductivity 
$*«0).  The  following  sequence  outlines  the  derivation  of  the  wove 
equations. 
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e 


b 

"  St* 

V  xK  =;  + 

meo&e'  +  &P' 

£  (?  *>?)  -  ?  x  (jj-j n 


“  r  x  v  xf 


1  V?  7*  2  -*\ 

=*  -  “  (v7  v  •  f  -  V  £  )  . 

Thus,  . 

dt2  St2 

Because  the  medium  Is  assumed  to  be  Isotropic,  propagation  can  be 
restricted  to  any  axis  without  lo68  of  generality.  Consequently, 
let  a  uniform  plane  wave  propagate  down  the  z-axls.  The  wave  equation 
has  now  become: 

d 2  02  ->  a2  —• 

p  £  ^eo — 2  ^  ^ — 2^  ' 

dz  °at  at2 

The  wave  equation  as  given  above  possesses  no  coupling  between  coordi¬ 
nates  due  to  vector  operations.  Therefore,  the  expression  of  the 
wave  equation  for  the  components  In  each  coordlantc  axis  Is  Independent 
of  the  other  tv/o.  It  Is  possible  to  examine  either  the  x-axls  or  y- 
axls  component  and  not  lose  any  accuracy.  Consequently,  no  distinc¬ 
tion  will  be  made  and  only  one  component  will  be  used  with  the  result 
that  the  wave  equation  Is  now  a  scaler  equation  of  the  form: 
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^£(z,t>^c  ~£(z,t)+u-^P(z,t)  . 

az  ar  atr 

Nov;  that  a  satisfactory  wave  equation  hao  been  obtained,  the 
relationship  between  the  exciting  field  and  the  generated  field  must 
be  obtained.  To  do  this,  let 
£(z,t)=£r(z,t)+£g(z,t)  . 

This  relates  the  instantaneous  field, £  (z,t),  to  the  fundamental  or 
exciting  field,  £^(z,t),  and  the  generated  field,  £g(z,t).  The  wave 
equation  now  becomes: 

o  2  2 

*“?  [£f(z,t)+6  (zft)}=ve  -~[£f(z,t)+£ _(z,t)}fji“^f>(z,t)  . 

3z^  1  8  °at  8  3t 

It  i 8  desirable  to  relate  these  functions  to  their  frequency  domain 
counterparts.  For  this  reason,  the  following  definitions,  previously 
preiented,  will  be  applied  to  the  latest  form  of  the  wave  equation. 


uu 

J  Ef(2,B)e-3k(»)zeMda, 

mCD 

GO 

EEU<tH?  f  E6(z.a>)e-Jk(<u)ze*«tdu) 

•«CO 

CD 


7 


Thus, 


*  00 


i'j  p(z,a.)etotcb]  . 


"erfcrming  the  indicated  differentiation  results  in 
<) 

2 

^  f  Ef(z<tu)“J2k(‘«)'^Ef(zitu)-k(iu)2E1.(z,a))]e“^k^^ze'^,t'da) 


2 

J  C_22Eg( z.«' ) - j2k('J) )^Eg( 2,® ) -k(io)  ^ z,«j)]e‘3k(“>)zei"tcb= 


•"otsj  U2)Ef(z,a)e-JkW2e>S^J 


00 

(nu^wJe^]  . 

-CD 

Assume  that  the  conversion  of  the  fundamental  to  the  third 
harmonic  is  sufficiently  small  that  the  amplitude  of  E^(z,«>)  is 
essentially  unchanged  as  the  wave  passes  through  the  medium.  Then 
the  derivatives  with  respect  to  z  of  E^(z,oi)  may  be  set  equal  to  zero. 
Using  the  above  assumption  and  dropping  the  common  factor  of  mj— ,  the 
wave  equation  is 
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-j  k(tu)^E^(  ZjuOe""^10^  Ze'J wtdurf-J*  [-^3g(  z,u))-j2k(ou)-j^g(  z,tu)- 

m£0  -£0 

CD 

-CO 

CD  CD 

-HC0j*  U)^Eg(  z,(o)e  3k(uj)  Ze^daj-v-J*  JJ^P(  . 

•40  -CD 

JL 

The  propagation  factor,  k(w),  is  equal  to  (nt)2ut  where  e*e(u>)  for 
most  materials  and/i=/uQ  (assuming  dielectric  materials).  Using  the 
relationship,  the  wave  equation  simplifies  to 

CD 

f  z>(u)-j2k(u))-|:E  ( z,u>) ]e~^k^ ze‘^tda) 

J  fa*  6  dz  g 


smVj  0)^P(z,u>)e^wtdu}  . 

-40 

This  may  be  rewritten  to  take  advantage  of  the  common  variable  of 
integration  and  an  equation  consisting  of  one  integral  can  be 
obtained. 

CD 

/2 

82  ,  , 

It  is  felt  that  the  factor  — rE  (z,w)  may  be  equated  to  zero  by 

dz  S 

assuming  that  the  envelope  of  the  wave  is  slowly  varying. 

« 

This  results  in  the  equation, 

CD 

f  [M.a)2P(z,u);-j2k(io)~E  . 

J  cz  g 


•40 
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For  this  relationship  to  be  true  over  the  nearly  unlimited  conditions 
placed  upon  it,  the  integrand  must  equal  zero.  Thus,  the  following 
equation  is  obtained 


Recall  that  k(u)=(/u)Iu  and  the  above  equation  becomes 


^g(z,U))=-j^^P(z,to) 


2e 


jk(ou)z 


At  this  point,  an  equation  relating  E  (z,io)  to  a  generating  function, 

6 

P(z,u)  has  been  obtained. 

C.  P(z,o)  )  os  a  function  of  Ej(z,m).  In  this  section,  a 

relationship  is  developed  between  P(z,u)  and  E^(z,«o).  In  IIA,  the 

following  definition  was  given: 

00 


<P(z,t> 


2n 


P(z,ou)e^d(ju 


J 

mJCO 


It  follows  that 

CO 


P(z.w)- J(P  (zJt)e""^dt 


Let(P(z,t)  be  defined  in  such  a  way  that  it  possesses  no  linear 
susceptibility  term.  This  may  be  achieved  by  allowing  the  permittivity, 
eQ*  of  the  proceeding  section  to  be  replaced  with  t  which  contains 
the  linear  susceptibility  term.  Because  the  interest  of  this  paper  is 
THG,  only  the  non-vanishing  fourth  rank  component  of  the  electric 
susceptibility  tensor  will  be  considered.  Thus, 

(P(z,t)=e0  Xff(z,t)'5 

where  X  is  the  non-vanishing  susceptibility  component. 


/P(z,t)=c0x[^y 


00 

5x(^)3C :fzf(z  jW '  )e"Jk(tu ' }  V“ '  **,  •  ]  . 


[ 


)e“Jk(u,")zeJ. 


00 

i-'"tcb"]C J 


E.(zjtu'")e";)k(tu",)zeJ''J,Mbc 


P(z,w)  Is  merely  the  Fourier  transform  of  this  equation. 
00  00 

P(z,whJ  {eQx(  ~)3J  Ef(z,U)')e“jk(uj')zeJ‘w't^'- 


00 

00 

*■00  «/„ 


This  may  be  rewritten  as 


P(^)=c0X(^)3 


Ef(z,(U1)Ef(zJa)")Ef(z,aj'")' 


■O  UO  UD  — ra 


e-jcwo-j+kU'o+kU"):  vj[«^  ■  vv]^  1<VW1W. 


The  relationship 
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may  be  used  to  obtain  the  equation 

03  CO  QO 

p(2,»K0x(i-)2_/’//Ef(  Z,» '  )Ef(  z,„")Ef(  z,w"%(ww  'VV'O- 

•^O  —m 

e-«^+k(»'O+k(»"0]^v  _ 

P(2,oj)  ■  0  except  when  '  r-.  0. 

This  condition  is^uced  in  the  form  ai»  to  get 

—CO 


e“j[k(conU,,-u),,,)+k(0'5xk((U'',)3zd;u,)dbl,) 

The  expression  derived  in  the  T"’-eceeding  paragraph  is  quite 
general.  At  this  point,  let  it  be  assumed  that  the  spectrum  of 
Ef(z,o.)  is  highly  concentrated  about  MqI  the  fundamental  frequency. 
The  degree  of  concentration  that  is  desired  is  enough  so  k(u)  may  be 
represented  by  the  first  two  terms  of  a  Taylor's  Series  Expansion. 
Under  these  conditions,  the  following  expression  is  valid. 


k(co)-k(u)  )+(m-u)  )~ 
°  O  0(1) 


Let  k(u)  )-k 
0  0 


U>=U), 


u^o 

Then 

k(co-u)  "■ -V " ")^ko-k-v  (to -tii 1 1  -tu '' Su  ) 
k(u) "  )=kQ-tar  (to  "-u)0) 
k((ju,n)=koio'((ju,"-»0)  , 

and  the  sum, 


and 

CO) 
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k(a)-tu%J"'i)+k(tu")+k((0l")=3k  Mw-yu  ) 

•ni _  0  0 


Thus, 


P(z,u>)-eox(^)y  J Ef(z>u)n0"-^,")Ef(z,u,")Ef(zJu,"')e'J  3Vjr(j,"3',5o)  ]zda)„dajn 


-OD  -CD 

Because  the  phase  factor  in  the  integrand  has  no  dependence  upon 
either  to"  or  o'",  it  can  be  removed  from  the  integral  and  placed  with 
the  constant  factor.  If 

CO  CO 


G (z,w)-J  J Ef(ziu;nol,^l")Ei:(zJa)")Ef(z>u),")djD"d;i;"'  , 

-CO  —CO 

then 


P(ZjU))=e 


1  2  -j[3k  )] 

e 


z 

G(z,tu) 


It  is  interesting  to  note  that  G(z,«)  is  the  Fourier  transform  of 

£f(z»t)  without  its  associated  phase  factor.  Because  6^(z,t)3  is 

essentially  unaltered  by  its  passage  through  the  non-linear  medium, 

3  3 

then  ^(ojt)  and  G(o,w)^G(z,w)#  It  con  also  be  seen  that 

G(o,u)  has  as  its  inverse  Courier  transform  £f(o,t)^.  Thus,  G(o,o>)<-> 
£f(o»t)  .  At  this  stage  in  the  derivations,  the  relationship  between 
the  fundamental  field  in  the  time  domain,  £f(z,t),  and  the  rote  of 
growth  of  the  third  harmonic  in  the  frequency  domain,  -~E  (z,«»),  is 

OZ  g 

emerging. 

D.  Determination  of  E  (z,m).  The  actual  spectral  density  of 

© 

the  third  harmonic  can  be  calculated  by  assuming  it  to  be  zero  at 

z=0  and  integrating  to  z. 
z 

Eg(z,(o)=J^Eg(z,tu)dz  . 


0 
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The  result  of  IIB  was  that 


Then 

z 

Eg M-^f PM°3kMziz  . 


Using  the  result  of  IIC, 
z 

e - 31^ J 


2  “JC3ko+o'(co— 3au^)  ]z 


)  e 


°  G(o,:»)}oJkWZd2  . 


Because  of  the  origin  of  P(z,w),  it  is  sharply  peaked  about  3a>o>  and 
the  expansion  of  k(w)  about  3w  is  valid. 


3k 


k(w)^k(3cu0)+U-3u)0)“ 


ur-3u)y 


Let 


dk 

3to 


(«= 3«>, 


then  k(ui)  =  k~  +  7  (<n-3<u  ) 
j  o 

and  then 


.  .  k(ou)e  x  r 

Ep,(z,cu)=  -j  2  I  G(o,(u)e 
g  2e(2n;  I 

o 

k(to)e. 


-j[3ko+of(o)-3a)0)]z  j[k^-+Y  (u>— 3u)q)  ]z 


dz 


=  0(o,„)  f 

2e(2rr)^ 


dz 


«  -j- 


k(w)e  j 

j - 2. 

2c(2n)J 


GCojtu) 


B~j[3k0-k^+(o'-y)(ai-3o>0)]z 


-j[3k0-k^+(a-y)  (o>-3oj0)  ] 


o 
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Evaluating  this  and  recognizing  that 


E„(z,U))  =  -j- 


2er(2n)' 


-J[3k0-ltj+(o/-*I')(aj-3^o)  ]z 
:°,yj)e-j[3k0-k3+(a-v)(u)-3u)0)  ]  * 


The  last  factor  can  be  converted  into  the  from  e  sin  Also, 
since  E^(z,«j)  is  highly  concentrated  about  3o>o,  let  the  amplitude 
factor,  k(w),  be  replaced  with  k^.  The  final  equation  for  this  section 


is  then 

E  (z,w)= 


2er(2n) 


2  G(o,uj)e 


-j[3k0-k3-h(o/-Y)  (u)-3ajQ)  ]z/2 


8in{  [3kQ-k3+(cy-Y)  (q)-3a)0)  ]z/2} 

[3k0-k3+(o'-Y  )  (aj-3cuQ)  ]z/2  * 

E.  Determination  of  F  (z,t).  The  final  step  is  to  convert 

ri.  r-  -  "  r  '  i  g 

E  (z,w)  into  its  time  domain  counterpart.  This  is  done  by  taking  the 
8 

inverse  Fourier  transform  which,  due  to  the  complexity  of  E  (z,w),  is 

8 

not  a  simple  task.  The  equation  which  gives  £  (z,t)  is: 

CO  8 


This,  upon  insertion  of  the  developed  form  of  E  (z,w),  becomes: 


2er(2n) 


Xz  -j[3ko-k34(a?-Y-)(u)“3uj  )]z/2 

- 5"  G(o,w)e  J 

A  \  ^ 


.  ^lC3l<0-V(o-Y)(i»-3a)0)>/2i  _jkW,  M 


[3kQ-k3+Ca'-Y  ) (urO^jz/  2 


eJ  d< 


The  integrand  contains  three  principal  factors  which  are  listed  here: 


l)  G(o,w), 


2)  e 


-j[3ko-k3+(a-\' )  (u)-3u)0)  >/2  sln{  [S^-k^+Cor-y)  (gj-3(o0)  ]z/2] 


[3k0-k3+(<y-v)(tD-3a)0)]z/2 


it  -  jk(iu)z 

,  assuming  that  k3,  x  and  ^  are  constant.  These  will  be 
treated  in  order. 

Earlier  it  was  stated  without  proof  that  G(o,c)  «_> £  (o,t)3.  Thi 
is  quite  true  except  for  the  scale  factor  associated  with  £  (o, l)3. 
The  derivation  of  this  relationship  follows, 

i  o  ”j[3k  +or(w-3'ju  )  ]z 
P(z,io)=eox(2^)  e  G(z,w) 


from  lie.  At  z=0 

P(o,m)=eox(^-)2  G(o,«) 

It  will  be  recalled  that 

(P(z,t)=eox£f(z,t)^ 

Therefore,  if2/(o,t)  «-*G<o,«),  then 

and  then,  by  simply  rewriting, 

^(o,t)=(2n)2£f(o,t)3  . 

The  second  of  the  three  principal  factors  is  the  most  difficult 
to  evaluate.  For  this  reason,  the  actual  development  has  been 
placed  in  Appendix  A.  The  result,  however,  is 


3k -k. 


'  jQf-VJ  Z 


,  o<t<(o'-v)z 


'  o,  all  other  t 

when  it  is  assumed  that  «-y>0.  A  change  in  the  sense  of  e-y  merely 
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results  in  a  change  in  the  sense  of  the  inequality. 

The  last  factor  is  Using  k(«.)=k,+y(„-3u,  ),  the 

3  o  ’ 

factor  becomes 

■"jr^v(a)-3o)0)]z 

which  is  equivalent  to 

■-jLko“'3u)  , 

e  ^  0  e-jyau 

The  last  form  consists  of  a  phase  factor  for  propagation  of 
e-J[l<3“3M0Y]z 

and  a  time  shift  of  7 z  seconds. 

The  final  bit  of  work  consists  of  using  the  information 

gathered  about  the  three  integrand  factors  to  obtain  £  (z,t).  This 

8 

hinges  on  the  fact  that 

F(u))G(cu)  *-^f(T)g(t-T)dr  . 

Identifying  F(w)  with  the  sin  <o/<d  function, 

G(«)  With  £f(o,t)3 

and  the  time  delay  with  the  dependent  variable,  the  following  ex¬ 
pression  can  be  obtained 


£e(s,t-^)3-  -  -2- 

6  2sr(2n)2 


,(a-Y)z  3k -k~ 

1  ^T-CV^oY)2) 

lor— yi  z 


•(2n)^££.(o,t-T)'3dT  . 

Define  a  change  of  variables 

T  '=t+yz 

and  t  '=^t+YX  • 

Then,  using  these  changes,  and  cleaning  up  some  simple  algebra, 
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JM 


otz 


V-’^  / 

Yz 


j£C3u,o“ 


£f(o,t '-T'^dT 


Dropping  all  of  the  primes  gives  the  third  harmonic 


yz 


,'i  I  e  £f(o,t-T)3dT  . 


In  the  time  domain,  the  privilege  of  examinl.:**  both  the  real  and  the 
imaginary  components  of  a  function  docs  not  exist.  Consequently,  a 
choice  must  be  made.  If  the  fundamental  function,  f^(o,t),  is  a  real 
function,  then  the  generated  function,  £  (z,t),  must  also  be  real  in 

s 

order  to  be  observed.  Applying  the  restriction  that  only  Re  {£  (z,t)| 

& 

i ®  meaningful,  then  £  (z,t)  for  this  paper  becomes 

8 


j(C3o)0. 


<*-Y 


]t+[ 


^2 

Qf”Y 


£^(o,t-T)^dx 


This  becomes  the  final  equation  of  section  II. 

otz 


-M 

£g(*'t)as32 


Yz 

^as  been  found  as  a  function  of  the  fundamental  field, 

u  1 

6f (o,t)  and  parameters  of  the  medium  such  as  k^»  '§»'*''***  z»  f r  ond  *• 


SECTION  III 


EXAMPLES 

A.  Eternal  Sinusoid.  To  illustrate  the  use  of  the  previously 

obtained  equation  and  as  a  basic  check  on  the  validity  of  the  theory, 

the  elementary  case  of  the  eternal  sinusoid  is  considered.  Let 

£,(o,t)*£  sin  u  t. 
f  ’  o  o 

Then, 

az 

x  f  3k -k_  3k  -k, 

EE(2'°-  ]  !il[(1'onM^]‘1 

yz 


•  £q2 3  sin3{u)o(t-i)}dT  . 

3 

Using  the  identity  sin  X=3/4  sinX-1/4  sin  3X  and  retaining  only  the 
third  harmonic  term  rctults  in 


Vz't} 


az 

r 


sln{  [3w0- 


3k  -  3k  — k_ 

]z}  sin{  3j»0(t-T)  }dT  . 


yz 


2 

When  the  group  velocities  arc  matched,  and  the  group  velocities 

equal  the  average  phase  velocity,^  )/ 2,  then 


£g(z,t)=- 


*3*e03* 

8c_ 


sin- 


cos(>,.)ot-k^z} 


2  .  Sec  Appendix  D. 

3  .  Sec  Appendix  C. 
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This  Is  Ir  excellent  agreement  with  previously  established  results^ 
which  have  the  Intensity  at  3u»o  proportional  to 


X 


2 


sin2(-Ak) 

1R 


Thus,  It  is  reasonable  to  assume  this  theoretical  development  to  be 
valid. 


B.  Gaussian  Pulse.  This  portion  of  Section  III  consists  of 
an  idealized  pulse  in  the  form  of  a  6ine  wave  which  is  amplitude 
modulated  with  a  gaussian.  The  third  harmonic  that  is  generated 
is  examined  in  both  the  time  and  frequency  domains  under  conditions 
of  both  matched  phase  velocity  and  matched  group  velocity. 

The  input  or  fundamental  pulse,  (o , t) ,  is  a  gaussian  modulated 
sine  wove.  The  temporal  characteristics  ore  completely  described  by 
the  equation 

2  2 

£f(o,t)«£oe"'t  /T  sin(<«ot) 


where  T  establishes  the  pulae  width.  The  spectrum,  E^(o,u>),  is  de 
fined  by  the  Fourier  transform  of  £^(o,t).  This  is  given  by 


Ef(o,tu)»‘3{£oe"t 


2/t2 


3in('juQt)} 


3{e-t2/T2}*3{sin  .joot} 


where,  os  usual,  3  |  |  denotes  the  Fourier  transform  and  the  astrisk 


4.  Bey,  Giuliani  and  Rabin,  "Phase-Hatched  Optical  Harmonic  Genera¬ 
tion  in  Liquid  Media  Employing  Aromalous  Dispersion",  IEEE  Journal 
of  Quantum  Electronics,  Vol.  QE-4,  number  11,  pp  932-939, 

November  1968. 


A 


(*)  denotes  convolution.  Appendix  E  is  a  proof  of  the  convolution 
relationship. 


V0<“>-tr  t/^T  e'T2"2A)*(jn[6(«A.0)-6(<.^.0):) 

,/7  f  Tre'T2(‘"~J’°)2A 

=j^2~  coLe  -e  ]  . 


The  power  spectrum,  S_(o,<..),  is  given  by 


S^,(c,tu)=Ef(o,a))Ef( 


°'t")  ^ 


where  the  asterisk  (*)  os  a  superscript  denotes  the  complex  conjugate. 
Thus, 

s(n  N  tt  _2  .  -t2(w-w0)2/2  -T2(u^  )2/2  2 

T  w  [e  +e  IE* 

if  it  is  assumed  that  the  positive  and  negative  frequency  components 
do  not  overlap.  This  assumption  causes  a  component  about  zero 
frequency  to  be  dropped,  but  in  the  case  of  a  loser,  this  causes  no 
difficulty. 

Next,  the  generated  field,  £  (z,t),  will  be  calculated.  From 

8 

Section  II. E,  the  following  equation  is  obtained: 

az 

lc,  X  r  qw-  -ib- 


E  (2,t)=^ — 2 - 

g  2  |or-Y  i  c 


Y* 

For  £f(o,t)=£Qe  t  /T  6in(m  t)  , 


•  sin3U  (t-r)]dT  . 


21 


*3 

By  use  of  the  trigonometric  identity,  sin  x-3/4  sin  x-1/4  sin  3x, 
in  the  above  expression  and  retaining  only  the  term  containing  3x 
because  the  only  interest  is  in  the  third  harmonic,  the  following 


equation  can  be  obtaine.d. 

„  ore 


Vz 

•  sinO^t-'OjdT  . 

The  identity,  sin  x  sin  y  =  [cos  (x-y)  -  cos  (x+y)]/2,  may  be  used 
to  rid  the  integrand  of  the  product  of  trigonometric  functions.  Thus, 


Vz,t)~l6 


lo-yi e 


-  oz 

J  e-3(t-r)2/T2 


3k  — k„  3k  -k 

•  c°sLfr’oi%r 2]T+['Pr2'  Y-k,]^->»0t)dT 
02 

-  [  «.-3(t-T)2/T2  .  rr  3kc-k3,  Jko-k1 

J  custt-  a-y  Y-kj]!H-3«)0t)dT  . 

yz 

For  the  case  where  the  phase  velocities  arc  matched,  e.g. 

3^0C^3  ^sce  Appendix  C  ),  the  expression  for  6  (z,t)  becomes 

8 


e  (,  t)  ye°  f  e-3(t-T)2/T2 

yz 

oz 

C  2  2 

•  cos[6JjoT-k3z-3o)0t}dT  -cos^t-k^z }/  e-3^""^  /T  dT 

rr  •  • 


For  the  case  in  vhich  the  gaussian  modulation  varies  slowly  v;hen 
compared  to  6wq  and  the  interval  (a-7)z»6wo,  the  first  integral 
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can  be  neglected  and  then 


az 


g^^^l6ia-y ,  c. 


C0s{3: 


e-3(t-r)2/T2  dT  _ 


\'z 


The  spectrum  fa’-  the  above  expression  is  given  by 

Eg(  *,«>)- 3  (£g(i>,t)) 


I  -jkjZ  -T2(ffl-3«i0)2/12  -j('»-3'JJ0)vz 


32/3  i  o-yi  e, 


oj-3'jj. 


+e 


jy  -T(m3u;0)2/12  e-^wf3o)0)°fZ  yJ(<”+3u>0)vz 
6  •  m+3a)„ 


The  pov.’er  spectrum  is  found  as  before.  Neglecting  the 
positive  and  negotive  frequency  overlap, 


S  (z,(o> 


k32X2fQ6T2nz2 

3072  e  2 


r  sin{(u)-3ooo)(ar-v)f} 
(ou-3tu0)(a-Y)| 


7°  "2  — T2  ( to  — 3'ju  Q) 2/ 6 
e 


T sin{  (m+3yu0 )  (a-\ )  |}" 

2  -T2(^0)2/6 

L  («ri3u)0)  (a-Y  )§ 

c 

The  above  calculations  are  similar  in  principle  to  those  performed 
earlier  for  the  fundamental  wave.  Due  to  the  more,  complicated  func¬ 
tions,  these  calculations  arc  more  involved  and  not  particularly 
enlightening.  Consequently,  only  the  results  arc  shown  here. 

The  third  harmonic,  £  ( z,t ),  can  also  be  found  under  conditions 

8 

of  matched  group  velocities  between  the  fundamental,  e^Cojt),  and  the 

third  harmonic.  The  output  under  matched  group  velocities,  c.g. 

a~y  (sec  Appendix  B),  is  obtained  by  taking  the  limit  of  £  (z,t)  os 

S 

a  approaches  7 .  The  treatment  is  as  follows. 
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6  (z,t)=liin 
5  O'-hY 


a  z 


l6\y-\ 


-3(t-x)2/T2 


[/ 


yz 

3k  —  k  3k  —  k 


a  z 


-  J  e-S^-)2/!2 


Yz 

Because  the  limit  Is  being  taken,  ic  is  reasonable  to  extract  the 
gaussian  function  from  the  integrand  by  substituting  7 z  for r  .  The 
equation  then  becomes 

z  V  ' 


e-3(t-v,)W  ^ 


r  a’z 


I vi 


f  3k  -It, 

j  cos{[6jj^-  ^]t 


O  Q'-'/ 


3k  -k 

+[— 7^  Y-k-jlz-Su^tjdT 

az 


yz 


r  3k  -k  3k  -k„ 

-  J  C0=ft'^p]’+[~ pV-^l^tJdT 


Yz 


This  can  be  reduced  to 

3k  -k- 

-k~x£  3z  sin{— f--^z}  3k  +lc, 

- i 2 —  - - -  ensf — ~ — 


Z-—  e*r  3VS 


£  (z,t)^ 


C  03  [  3'0ftt  -3to0Y  2  ^  g  -3  ( t  -Y  z)  2/T  2n .  m  o-v 


lim 

o— Y 


O'-*/ 


The  limit  in  the  above  expression  merely  establishes  the  sign  of 
£  (z,t).  A  sin  (x)/x  relationship  with  respect  to  the  phase  velocity 

O 

match  is  obtained  while  an  oscillating  relationship  is  obtained 
between  the  average  of  the  phase  velocities  and  the  group  velocity 
which  is  then  weighted  by  z. 
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The  various  spectra  are  found  to  be 
3  _  3k  -k. 

xJo2^El 

16/r  cp  3vki 

I  ■  — *  <J 


s 


.cost^ 


f  e 


(urf3co  )  / 12 


J  O^Y  a“v 


and 


lb 


v 


c-fv 


c03={^z.3,oYZ)je-T‘(->«o)^+e-TVj».)2/6j 


A  concluding  comment  and  calculation  give  some  insight  as  to  the 
effects  of  the  phase  velocities  and  group  velocities  on  TUG.  If  the 
peak  intensity  of  the  THG  pulse  is  calculated  and  then  normalized  to 
the  peak  intensity  obtained  when  the  independent  variable  is  zero, 
then  a  plot  of  this  ratio  can  be  made  and  a  qualitative  feeling  for  ’ 
the  process  developed.  In  the  case  of  matched  phase  velocity,  the 


ratio  is  given  by  the  equation 


I 


25 


/\ 


where  I  h  peak  intensity  as  a  function  of  the  independent  variable, 

peak  intensity  under  the  condition  of  matched  group  velocity 
(a =7)  and 

T  5  a  parameter  directly  related  to  pulscwidth. 

In  the  case  o£  matched  group  velocities,  the  expression  is 


?  Sln  *•  2 


3k«-k, 


i"  f3ko_k3  ,2 

m  I — 5 — *z] 


when  the  group  velocitier.  equal  the  average  phase  velocity.  These 
functions  arc  plotted  and  presented  in  the  accompanying  Figures  I 
and  II. 


C.  Linearly  Chirped  Rectangular  Pulse.  This  portion  of  Section 
III  presents  the  types  of  calculations  made  and  the  results  obtained 
from  these  calculations  for  the  third  harmonic  generated  by  a  linearly 
chirped  rectangular  pulse.  The  calculations  will  be  only  briefly 
described  because  of  the  complexity  of  the  expressions  and  the 
difficulty  and  bulkincss  involved  with  executing  the  required 
operations. 

Assume  a  fundamental  pulse 

eo cosuot+|^t2},  - 

'  o,  otherwise 

The  spectrum  of  this  pulse  is  given  by 


Ef(o,u))-9  [£f(o,t)}  . 


dt 


By  the  use  of  the  Fresnel  Integrals, 


Following  the  previous  examples,  it  is  seen  that  the  spectral  power 
density,  S^(o,cj),  is  equal  to 


The  equations  presented  in  the  proceeding  paragraph  arc  fairly 
large  and  do  not  leave  the  reader  with  a  reasonable  feel  for  the 
shape  of  the  function.  The  equations  developed  for  the  THG  as  a 
result  of  this  fundamental  excitation  become  much  moru  cumbersome  and 
unwieldy.  For  these  reasons,  information  on  the  spectra  will  not  be 
presented. 

The  third  harmonic,  £  (z,t),  that  is  generated  by  the  chirped 
fundamental  will  be  calculated  for  the  case  of  matched  phase  velocities 
between  the  fundamental  and  the  third  harmonic.  From  Section  II, 
port  E, 

-  5]^ir  ^ 

J  yz 

•  f^(o,t-T)^  dT  . 

Lnder  the  condition  of  matched  phase  velocities,  3kQ=k^ ,  then 


-k-X 

VZ,t)  =  2  la-vf 


Sotz 

sin[3'JU0T-k^z)ff(o,t-T)^  dT  . 
yz 


Because  5^Co,t)  is  expressed  in  a  piece-wise  fashion,  £^(z,t)  is 

also  obtained  in  o  piece-wise  fashion.  The  limits  of  the  integral  in 

£  (z,t)  are  affected  and  the  expression  of  £  (z,t)  is  rather  cumbersome. 
8  8 

An  explanation  of  the  effect  that  £^(o,t)  has  on  £^(z,t)  with  respect 

to  the  integration,  and  an  example  establishing  the  integral  limits 

for  the  various  regions  of  8  (z,t)  is  presented  in  Appendix  F. 

8 

There  is  a  general  form  for  £  (z,t),  under  various  conditions  of 

S 

a,  7,  z,  T  and  t.  The  following  equation  presents  this  form  while  the 
accompanying  table,  Table  I,  presents  the  conditions,  regions,  and 
integral  limits  for  the  four  coses  considered. 

In  general, 

-k f  2 

£g(z,t)  =  Qjalyi  °  J  sin{3u)oT-kQz}  cos[3cuQ(t-T)  +  |  (t-T)2}  d7  . 

T1 

E  (z,t)  is  necessarily  zero  except  in  regions  I,  II  and  III  as  given 
8 

in  Table  I .  The  actual  solution  for  the  first  two  cases  will  be 
given.  Thi3  is  considered  to  be  sufficient,  because  these  consider 
the  situations  in  which  the  pulse  is  smaller  than|a.-y|z  and  also 
larger  than  la-7lz.  Interchanging  the  group  velocities  does  not 
cause  a  drastic  change  in  the  character  of  the  solution.  The  solutions 


' 

I 

I 


arc: 


X^«-Y|cr  2]  •°[V?  ]sin  (3»0t^3*) 

+  *]  <t"/Z)]  ]C0S  (3*Ot'lt3Z)]  •  V2-|lt<Y2+| 


For  all  t  other  than  specified  above,  £  (z,t)  is  necessarily  zero. 


e-101 
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The  third  harmonic,  £  (z,t),  that  1*  generated  under  matched 

8 

group  velocity  conditions  will  now  be  calculated.  Thij  a  obtained 

by  evaluating  £  (z,t)  In  the  limit  as  a-*?.  Thus, 

8 

EgU,t)  1±.  21*!,.,  (  ain(O0  - 


•  £f(o,t-t)3  dT  . 


For  the  chirped  pulse  case,  thls-can  be  seen  to  reduce  to 


-lc,xf 


£g(z,t)~  —  costeju^t-Yz)  +  |  IT  (t-Yz)2}  lim 


jyz 

\  t+c 

*  \jrr 


3Vk 


Y-^>} 


Upon  evaluation,  £  (z,t)  becomes 
& 


3Vk3 

-l^xe^z  8int  2 

&r  “5*^ 

— 5 — 


(s-) 


3k  +; 

*1h{3cu0yz  -  —y 


}  cos[3o)n(t-Yz)+|  p  (t-vz)2] 


*g(z,t) 


T  T 

for  yz  ~2-^-^z  +  2 


0,  otherwise 
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The  effects  of  several  parameters  on  peak  Intensity  will  be 
shown  now.  As  before,  the  Information  will  be  shown  graphically  and 
In  a  normalized  manner.  The  development  used  Is  quite  similar  to  that 
of  the  gaussian  case,  although  this  time  It  Is  not  possible  to  present 

i 

as  much  information  as  would  be  desired  due  to  the  number  of  variables 
Involved  and  the  complexity  of  the  function.  For  the  phase  velocity 
matched  case,  the  particular  case  of  a>y,  T  <l«-riz  will  be  considered. 
In  addition,  only  Region  II  of  the  solution  will  be  evaluated  under 
the  assumption  that  It  domlnstes  the  peak  Intensity.  After  normalizing, 
the  peak  Intensity  Is  given  by: 

^C(X)2  +  S(X)2^/tt?,  where  X  . §  , 

thus  containing  both  n  and  T.  A  plot  of  this  function  Is  shown  In 

Figure  III,  where  the  pulse  width,  T,  Is  held  constant. 

Another  case,  when  JT  Is  held  constant  and  T  is  the  variable 

parameter,  is  also  presented.  The  peak  intensity  is  proportional  to 

2  2 

C(x)  +  S(X)  and  a  plot  of  this  function  is  shown  in  Figure  IV. 

Finally,  the  case  of  matched  group  velocity  is  show;  in  Figure  V. 
This  is  simply  the  sin  (x)/xZ  function  that  appears  to  be  the  charac¬ 
teristic  function  associated  with  mismatched  phase  velocities. 


SECTION  IV 
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CONCLUSIONS 

A  theory  for  THG  has  been  developed  from  Maxwell's  Equations  for 
a  uniform  plane  wave  in  an  isotropic  medium.  This  theory  is  fairly 
general  with  the  only  outstanding  assumptions  being  that  the  generation 
process  is  assumed  to  be  relatively  inefficient  and  that  the  spectral 
bandwidth  of  the  pulse  about  the  center  frequency  is  sufficiently  small 
that  the  propagation  coefficients  of  the  fundamental  and  third  har¬ 
monic  can  be  represented  by  a  truncated  Taylor's  expansion.  For  the 
present  state-of-the-art  in  optics,  these  assumptions  are  quite  prac¬ 
tical. 

A  series  of  three  examples  were  then  examined.  The  first  of  these 
an  eternal  sinusoid,  was  examined  to  determine  whether  the  results 
of  this  theory  fit  previously  obtained  results.  A  comparison  was 
made  with  the  results  of  a  recognized  group  of  workers  in  the  field 

and  a  favorable  outcome  was  obtained. 

The  second  example  was  a  sine  wave  that  was  amplitude  modulated 
with  a  gaussian  wave  form.  This  fundamental  pulse  was  then  examined 
for  THG  under  conditions  of  matched  group  velocities  and  of  matched 
phase  velocities.  For  matched  group  velocities,  the  variation  of 
peak  intensity  followed  a  sin2  (x)/x2  function  while  for  matched 
phase  velocities,  the  peak  intensity  followed  a  erf  (x)2/x2  function. 

The  third  example  was  a  linearly  chirped  rectangular  pulse,  itiis 
example  was  quite  difficult  from  a  mathematical  point  of  view  and 
only  the  most  interesting  results  and  calculations  were  presented. 

The  THG  pulse  was  examined  under  the  conditions  of  matched  gr°uP 
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velocities  and  of  matched  phase  velocities.  For  the  condition  of 

matched  group  velocities,  the  peak  intensity  followed  a  sin2  (x)/x2 

function.  For  the  condition  of  matched  phase  velocities,  the  THG  pulse 

was  examined  as  a  function  of  frequency  modulation  and  a3  a  function 

of  pulse  width.  In  both  cases  oscillatory  behavior  was  observed. 

When  the  frequency  modulated  case  was  examined,  the  peak  intensity 

followed  a  (C(x)  +  S(x)  )/x  function.  When  the  pulse  width 

modulated  case  was  examined,  the  peak  intensity  followed  a  C(x)2  + 

2 

S(x)  function.  This  oscillatory  behavior  is  quite  interesting  and 
could  possibly  be  seen  in  an  experiment  if,  for  example,  a  train  of 
pico-second  laser  pulses  were  to  change  pulse  width  as  the  pulse 
train  progressed  in  time. 

In  conclusion,  this  work  has  developed  a  seemingly  valid 
theory  of  THG  and  has,  within  limits,  begun  to  identify 
THG  results  as  a  function  of  material  and  excitation  parameters.  It 
is  sincerely  hoped  that  this  effort  will  benefit  others  in  this 
field  of  endeavor. 
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Case 

Region 

Interval 

'1 

'2 

a  >  v 

I 

T  T 

yz  -  ^  t  <  az  -  ^ 

Yz 

T 

t  +  2 

T  >  |a-Y  |  z 

n 

T  T 

az-^<t<YZ+2 

YZ 

az 

m 

T  T 

Yz  +  2-  t-az  +  2 

1  2 

az 

a  >  Y 

i 

T  T 

Yz-£<t<Yz  +  2 

Yz 

t  +  I 

T  <  |a-Yl  z 

n 

T  T 

Y  z  +  ^<t<az«2 

t  .I 
t  2 

t  +  f 

m 

T  T 

az-2<t<az  +  2 

*-I 

az 

Y  >  a 

i 

T  T 

az  -  -  <  t  <  yz  -  2 

t+| 

az 

T  >  |  cr— V  |  z 

n 

T  T 

YZ--<t<QfZ  +  - 

Yz 

az 

in 

T  T 

QfZ  +  ^<t<Yz  +  2 

Yz 

t  2 

Y  >  a 

i 

T  T 

az-  —  <t<o'Z  +  — 

T 

t  +  | 

az 

T  <  |a-Y|  z 

ii 

T  T 

az  +  -  ^  t  <  yz  -  g 

T 

t  +  2 

t  .1 
t  2 

m 

T  T 

YZ  -  g  <  t  <  yz  +  j 

Yz 

t  .1 
2 

TABLE  I.  Integral  Limits  for  THG  of  Linearly 
Chirped  Pulse. 
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EVALUATION  OF  THE  INVERSE  FOURIER  TRANSFORM  OF: 
-J[3k0-*3+(«-v)  (oj-3o)0)  ]z/2  sin{  [3kQ-k3+(of-v)  (cu-3o)Q)  ]z/2} 
1  [3k0"k3 4^0f"^  )  Cod— 3cu0)  ]z/ 2  ~ 


Postulate  a  function  f(t)  and  knowledge  of  its  Fourier  transform, 
F(w).  Denote  the  transform  pair  relationship  by  the  symbology 
f(t)  < — »F(w) 

Let  F((u)  .  storA^-^l 

A(u)-yl) 

From  the  expression  In  the  title,  it  is  seen  that 
A  s  (cHY)z/2 


*0*1 


a  =  3iu.  - 

o  a-Y 


Thus  f(t) 


e-JA(u jra)  sinfAQo-.^)) 

A((o-a] 


The  scaling  relationship, 

f(rt)  —  Jfj  F(iS)  , 

is  applied  yielding 
f(At)  =  -i-  ej(tu"A)  szsbtal 

IAI  0U“A 

The  frequency  shifting  property. 


F(cu-u)0) 


f(t)  e 
is  applied  yielding 

f(At)  e_inAt  ♦->  -i-  e“^  sjn  .j, 
IAI  (0 

Finally,  the  time  shifting  property, 


f(t-t  ) 

o 


F(w)  e"Jwt° 


is  applied  and  the  scale  factor  IAI  is  transposed. 


Thus, 
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If  ,  then 

d(t)  -  (%,  -1  <  t  <  1 

<  0,  otherwise  . 

There  now  exists  the  equivalence 
d(t)  -  !AI f (A[ t+l])  e"jnA(t+l)  . 

Allow  a  change  of  variable, 
t»  -  Af t+l]  . 

Thus  f(t»)  -  d(—y  -  1)  e*"*'  . 

Deleting  the  primes  and  using  the  definition  for  d(t), 

4(0  -  jw e)Qt  • 0  £  ‘  i  »  ■ 

*0,  otherwise  . 

If  A  <  0,  the  inequality  is  reversed.  Using  the  definitions  of  A  and  n 
originally  stated, 

/ 

V  3k  -k, 

1  J[  3»  -  °  ]  t 

7^7  •  °  "r 

f(t)  - 

0,  otherwise  . 
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EVALUATING  £  (z,t)  WHEN  <*-Y 

O 


Let  £^(o,t)  =  sin  (DQt 

r  ,  .  "V  fott  3k  -k,  3k  -k- 


•  £o^8ln^t“»o(t-T)3dT  . 


Considering  only  the  third  harmonic,  £  (z,t)  becomes 
C«(,'°  "  0lnt[31"0'  ^]T+C^? 


•  sln{3u)o(t-T)}dT  . 


Using  sin  x  sin  y  a  |[cos(x-y)-cos(x+y)], 

£g(z't)  *’  li.k,-,°i€r  costc^-  -sh.y^-^1  Y-kjJz-^tJdT 

-/: C03[[-  -^:^3»0t)dT]  . 


Performing  the  indicated  integration  and  simplifying. 


U5 


laf  3  / 

fg(z,t)  *  16itt-vi  er  JiT^kT  I  »m(&»0M-31Co*->»0t) 

v-^r  \ 


-  aint^Yz-k^z-^t} 


3ko-k 


sin{  3a)Qt-3koz }  -sin[  ^t-k^z } 


Taking  the  limit  as  a-*y, 


v  >r  3 

£E(Zlt)  *  ^7"  3iv^  >in( <o>0V*-kjZ-3«)0t } 

-  «ln(fe0Ys-3lc0s-3ii0t)+sint3o!0t-3k0z]-slji(3(i)0t-li3z]j  . 


Using  sin  x  -  sin  y  «  2  sin  £(x-y)  cos  £(x+y). 


C  /  t,  V^c3*  ^HP*f 
V1-0  -  —  < 

—o  *  L 


cos(3<u0t+  —°2  J  -6u)0Vz) 


+cos(3u) 


*-3£ 


O  2 


z]  . 


Finally,  using  cos  x  +  cos  y  »  2  cos  £(x+y)  cos  £(x-y), 

,  V^o3*  z3  .  3kn-k, 

CgVZ,t)  «  “  3ko-ic,  cos(3uj0Yz-  — 2^  z}cos{3a)0t-3a)nYz} 


o  o' 
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REIATING  k  TO  PHASE  AND  GROUP  VELOCITIES 

The  accepted  definition  of  phase  velocity  ia  the  velocity  at  which 

a  fixed  poir, t  in  phase  travels.  For  a  simple  traveling  sine  wave, 

sin  (wt-kz),  the  phase  velocity  is  given  by 
u 

P  "  k 

Thus,  k  ■  —  . 

P 

If  3k  ■  k,,  then  -  ■  — 

03  ».  • 

and  since  -  3»  ,  the  conclusion  is  that  V  ■  V 

po  p3 

Historically,  the  phase  velocity  has  been  considered  to  be  the  speed 

of  light  in  the  medium  under  consideration.  Thus, 

VI  ■  c/n 
P 

where  c  »  speed  of  light 

n  -  index  of  refraction. 

If  Vp  "  Vpj »  t*ien  necessarily  n^  ■  n^  . 

Packets  of  electromagnetic  energy  are  usually  considered  to 
travel  at  a  velocity  called  group  velocity  rather  than  the  phase 
velocity.  Using  a  spectrum  closely  grouped  about  the  group 
velocity  la 


Thus  a  and  7  as  used  in  the  paper  are  reciprocal  group  velocities. 
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APPENDIX  D 

RELATING  TO  n  AND  & 
o«*  oA 

Bccfluae  moat  Information  relating  propagation  to  materiala  uaed 
in  optica  ia  given  in  the  form  of  n(A),  it  ia  appropriate  to  relate 
the  group  velocity  to  n(\). 


but  A  ■  AQ/n,  where  aq  ia  in  vacuo,  then 
Becauae 

2*  m 

ft  .  •  on 
“c  +  c*7  ‘ 

The  more  common  optical  term  la  diapcraion  which  la  Thia  will 

in  ®  o 

be  related  to 
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A  CONVOLUTION  RELATIONSHIP 


Prove : 


F(«)*G(«) 


Proof: 


-gr  F((fl)*G(«)  < 


f(t)g(t) 


^iM 


•u)G(u)du 


F(<u-u)G(u)du  e  d« 


^lo<u>  r(-“ 


)  e^wt  du  >  du 


Due  to  the  frequency  shifting  property, 

J® 

F(w-u)  e^**  dw  »  f  (t)  e^ut  . 


— ® 
Thus 


f  °° 

^  F((u)*G(ou)  «-*  -1-  I  G(u)f(t)ejut 

Ja> 


00 


du 


►f(t)  ^  •£—  1  G(u)  e^Ut  duj 
—00 


■£-  F(ou)*G(oj)  < - ►  f(t)g<t)  . 
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APPENDIX  F 

ESTABLISHING  THE  LIMITS  OF  INTEGRATION  OF  £  (z,t)  FOR  A 

8 

RECTANGULARLY  PULSED  £f(o,t) 

For  expressions  of  f  j(ojt)  that  are  in  a  piecewise  manner 
rather  than  a  single  analytically  expressed  function,  establishing 
the  limits  of  integration  for  E  (z,t)  is  not  alwaya  simple.  It  is 
the  purpose  of  thia  appendix  to  illustrate  a  useful  graphical 
technique  for  determining  these  limits. 

Let  £  (o,t)  -  (1,  - 


r* 

(0, 


2  —  2 
otherwise 


and  let 

az 

£,<*.<> 

7z 

Such  a  function  is  sufficiently  close  to  the  actual  function  of 

Cj(o,t)  and  £g(z,t)  for  this  illustration  to  be  valid  but  simple 

enough  to  convey  easily  the  intent  of  the  illustration.  Assume 

azyrz  and  that  T<£la-7|z.  For  t^yz-'j,  Figure  A-IV(a)  represents  the 

region  where  the  integrand  is  non- zero  and  the  region  of  integration 

do  not  overlap.  The  integral  is  necessarily  zero  for  this 

condition.  As  t  increases,  £^(o,t-r)  moves  to  the  right  and  finally' 

T 

enters  the  region  of  integration.  This  occurs  at  t»Yz-^  and  the 

T  T 

limits  of  integration  become  yz  and  t+j  until  t+^az.  This  condition 
is  shown  in  Figure  A-IV(b).  As  £f(o,t-r)  continues  to  move  to  the 
right,  finally  it  spans  the  entire  region  of  integration  and  the  inte- 

m 

gral  limits  become  yz  and  «z.  This  condition  begins  at  t-^oz  and 
T 

lasts  until  t-jp«}z.  Thia  condition  is  shown  in  Figure  A-IV(c). 


51 


After  t-^yz,  6f(o,t-r)  begins  emerging  from  the  region  of  integration 
T  T 

and  until  t-^a«z,  the  integral  limits  are  t-^  and  orz.  This  is  shown  in 

Figure  A-IV(d) .  The  regions  depicted  in  Figures  A-IV(b),  (c),  and  (d)  are 
called  Regions  I,  II  and  III  respectively  and  are  the  regions  in 
which  £  (z,t)  is  not  necessarily  zero.  For  t-^->az,  £..(o,t-r)  has 

O  A  I 

completely  emerged  from  the  region  of  integration  and  £  (z,t)  is 

g 

zero.  This  is  shown  in  Figure  A-IV(e). 

The  results  of  this  exercise  and  the  remaining  cases  involving 
«z,  yz  and  T  are  tabulated  in  Table  A- IV  for  an  integral  as  shown 
here: 

t2 

£g(zft)  -f  £f(o,t-r)3  dr  . 

Jti 
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6g(z,t)  -  J  £f( 0,t-T)3  dT 


-gas? 

Region 

Interval 

i 

12 

a  >  Y 

I 

T  T 

Yz-^<t<az-^ 

Yz 

t+l 

T  >  |ff-Y|  z 

II 

T  T 

az--<t<Yz  +  | 

Yz 

arz 

m 

T  T 

Yz  +  2-  t-az+2 

“I 

QfZ 

or  >  Y 

i 

Yz-|<t<Yz  +  | 

Yz 

T  <  |cr-V  |  z 

n 

T  T 

Yz  +  “<t<c*z-| 

“I 

m 

T  T 

0fZ  -  2-  t^az  +  2 

*-§ 

erz 

Y  >  a 

i 

az-|<t<Yz-| 

t  +  i 

arz 

T  >  |o-y  |  z 

n 

Yz-|<t<ffz  +  | 

Yz 

arz 

hi 

oz  +  ^<t<Yz  +  “ 

Yz 

1  “2 

y  >  a 

i 

az-|<t<crz+| 

t  +  i 

arz 

T  <  |  or— v  1  z 

n 

T  T 

orz  +  -  <  t  <  Yr  -  j 

t  +  i 

*-§ 

HI 

T  T 

Yz--<t<yz  +  | 

Yz 

*-a 

Table  A-I7. 

Convolution  Integral  Limits 


Mcv  EhIcv  E-*|CV  E~*|CM  .EhICM 


Eh|C\J 
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